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Abstract 

The nuclear matrix elements of the neutrinoless double beta decay (On/3/3) of most nuclei 
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is used to adjust the most relevant parameter, the strength of the particle-particle interaction. 
New results confirm that with such procedure the values become essentially independent 
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by reducing the coupling constant qa , as well as to the uncertainties of parameters describing the 
short range nucleon correlations. Theoretical arguments in favor of the adopted way of determining 
the interaction parameters are presented. Furthermore, a discussion of other implicit and explicit 
parameters, inherent to the QRPA method, is presented. Comparison is made of the ways these 
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I. INTRODUCTION 


Inspired by the spectacular discoveries of oscillations of atmospheric [1], solar [2, 3, 4, 5] 
, and reactor neutrinos[6] (for recent reviews see [7, 8, 9, 10, 11]) the physics community 
worldwide is embarking on the next challenging problem, hnding whether neutrinos are 
indeed Majorana particles as many particle physics models suggest. Study of the neutrinoless 
double beta decay {Qvj3j3) is the best potential source of information about the Majorana 
nature of the neutrinos [12, 13, 14, 15]. Moreover, the rate of the Oz//?/? decay, or limits 
on its value, can be used to constrain the neutrino mass pattern and the absolute neutrino 
mass scale, i.e., information not available by the study of neutrino oscillations. (The goals, 
and possible future directions of the field are described, e.g., in the recent study [16]. The 
issues particularly relevant for the program of 0z//5/3 decay search are discussed in [17].) 

The observation of Qvj3j3 decay would immediately tell us that neutrinos are massive 
Majorana particles. But without accurate calculations of the nuclear matrix elements it will 
be difficult to reach quantitative conclusions about the absolute neutrino masses and mass 
hierarchies and conhdently plan new experiments. Despite years of effort there is at present 
a lack of consensus among nuclear theorists how to correctly calculate the nuclear matrix 
elements, and how to estimate their uncertainty (see e.g. [15, 18]). Since an overwhelming 
majority of published calculations is based on the Quasiparticle Random Phase Approxi¬ 
mation (QRPA) and its modihcations, it is worthwhile to try to see what causes the sizable 
spread of the calculated values. Does it reflect some fundamental uncertainty, or is 
it mostly related to different choices of various adjustable parameters? If the latter is true 
(and we believe it is) can one hnd and justify an optimal choice that largely removes such 
unphysical dependence? 

In the previous paper [19] we have shown that by adjusting the most important parameter, 
the strength of the isoscalar particle-particle force so that the known rate of the 2z//?/?-decay 
is correctly reproduced, the dependence of the calculated Oz//?/? nuclear matrix elements 
on other things that are not a priori hxed, is essentially removed. In particular, we have 
shown that this is so as far the number of single particle states included is concerned, and the 
choice of the different realistic representations of the nucleon G-matrix. In [19] we applied 
this procedure to the Oz//?/? decay candidate nuclei, ^®Ge, ^°°Mo, ^^°Te, and ^^®Xe. 

In the present work we wish to expand and better justify the ideas presented in [19]. 
First, the method is systematically applied to calculate the nuclear matrix elements for 
most of the nuclei with known experimental 2z//?/?-decay rates. Second, the sensitivity of the 
results to variation of other model parameters is tested. These are the axial vector quenching 
factor, commonly described as a modihcation of the constant qa, and the parameters that 
describe the effect of the short range correlations. Finally, arguments in favor of the chosen 
calculation method are presented and discussed. 


II. DETAILS OF THE CALCULATION OF 0v(3(3 DECAY MATRIX ELEMENTS 

Provided that a virtual light Majorana neutrino with the effective mass {mpp), 
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is exchanged between the nucleons [20] the half-life of the 0z//3/3 decay is given by 

-L 1/2 

where G^‘'{Eq, Z) is the precisely calculable phase-space factor, and is the correspond¬ 
ing nuclear matrix element. Thus, obviously, any uncertainty in M'^'^ makes the value of 
{rriisf}) equally uncertain. 

The elements of the mixing matrix and the mass-squared differences Am^ can be 

determined in oscillation experiments. If the existence of the Qvj3j3 decay is proved and 
the value of T 1/2 is found, combining the knowledge of the hrst row of the neutrino mixing 
matrix \Uei\^ and the mass-squared differences Am^, a relatively narrow range of absolute 
neutrino mass scale can be determined, independently of the Majorana phases a* in most 
situations [14, 22]. However, such an important insight is possible only if the nuclear matrix 
elements are accurately known. 

The nuclear matrix element M'^'^ is dehned as 



where |i), (|/)) are the wave functions of the ground states of the initial (hnal) nuclei. We 
note that for qa = 1.25 nuclear matrix element M'^'^ coincides with of our previous 
work [19]. We have chosen this, somewhat awkward, parameterization so that we could later 
modify the value of qa and still use the same phase space factor G^'^{Eq, Z) that contains 
g\ = (1.25)^, tabulated e.g. in Ref. [23]. 

The explicit forms of the operators and are given in Ref. [23]. In order 

to explain the notation used below we summarize here only the most relevant formulae. We 
begin with the effective transition operator in the momentum representation 
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For simplicity we do not explicitly indicate here the q^ dependence of the form factor gA{q^) 
(the usual dipole form is used for both gv{,q^) and gA{,q^))i and the terms containing l/m^. 
The full expressions can be found in Ref. [23] and are used in the numerical calculations. 
Note that the space part of the momentum transfer four-vector is used in Eq. (5) since the 
time component is much smaller than q=\q\. 

The parts containing ^ + ml come from the induced pseudoscalar form factor gp for 
which the partially conserved axial-vector current hypothesis (PCAC) has been used. In 
comparison with most of previous Oz//?/? decay studies [24, 25, 26, 27] the higher order 
terms of the nucleon current (in particular the induced pseudoscalar as well as the weak 
magnetism) are included in the present calculation, resulting in a reduction of the nuclear 
matrix element by about 30% [23]. In the numerical calculation of M'^'^ here the summation 
over the virtual states in the intermediate nucleus is explicitly performed. 
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The relatively important role of the induced nucleon currents deserves a comment. In the 
charged current neutrino induced reactions the parts of the cross section containing gp are 
usually unimportant because they appear proportional to the outgoing charged lepton mass. 
That comes about because the equation of motion for the outgoing on-mass-shell lepton is 
used. Here, this cannot be done since the neutrino is virtual, and highly off-mass-shell. 

In order to calculate the matrix element in the coordinate space, one has to evaluate 
hrst the ‘neutrino potentials’ that in our case explicitly depend on the energy E™' of the 
intermediate state. 


/fx(ri 2 ) 


-ngX Jo 
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hK{q^)qdq 
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where K stands for E,GT,T. fp^Griqru) = joiqru) and friqru) = -j 2 {qri 2 ), where jo ,2 
are the spherical Bessel functions. The form factor combinations hx are dehned in Eq.(5), 
ri 2 is the distance between the nucleons, R is the nuclear radius, and Ei{Ef) are the ground 
state energies of the initial (final) nuclei. We note that in Ref. [23] the tensor potential is 
presented incorrectly with the spherical Bessel function jo{x) = sin(a;)/a; instead of j 2 {x). 
However, the numerical results were obtained with the correct expression. 

The individual parts of the matrix element , Eq.(3), are given by the expression: 


Mx= T. {-iy-+E’+j+^^2j +11 i.] 

J^,ki,kf,J pnp'n' [ dp' dd ] 

(p(l),p'(2); II f{ri 2 )OKf{n 2 ) II n{l),n'{2)-,J) 
x(0;i|[c+i.,],||J-A;/)(J-A;^|J-A;,)(J-% . (7) 

Here we dehne only those symbols that are needed further, for full explanation see [23]. 
The summation over J’’’ represents the summation over the different multipolarities in the 
virtual intermediate odd-odd nucleus; such states are labeled by the indices ki,kf when they 
are reached by the corresponding one-body operators acting on the initial or final nucleus, 
respectively. The overlap factor {J'"kf\J'"ki) accounts for the difference between them. 

The operators Ok,K = E,GT,T contain the corresponding neutrino potentials and the 
relevant spin and isospin operators. Short range correlation of the two initial neutrons and 
two hnal protons are described by the Jastrow-like function /(ri 2 ), 

firu) = 1 - - 72ri2) , (8) 

where the usual choice is [28] 71 = 1.1 fm^, 72 = 0.68 fm^. (These two parameters are 
correlated.) 

Finally, the reduced matrix elements of the one-body operators c+Cn (c^ denotes the 
time-reversed state) depend on the BCS coefficients Ui,Vj and on the QRPA vectors X,Y 
[23]. The difference between QRPA and RQRPA resides in the way these reduced matrix 
elements are calculated. 

As in [19], the quasiparticle random phase approximation (QRPA) and its modihcation, 
the renormalized QRPA (RQRPA), are used to describe the structure of the intermediate 
nuclear states virtually excited in the double beta decay. We stress that in the QRPA 
and RQRPA one can include essentially unlimited set of single-particle states, labeled by 
p, n in Eq.(7), but only a limited subset of conhgurations (iterations of the particle-hole, 
respectively two-quasiparticle conhgurations), in contrast to the nuclear shell model where 


4 



the opposite is true. On the other hand, within the QRPA there is no obvious procedure 
that determines how many single particle states one should include. Hence, various authors 
choose this crucial number ad hoc, basically for reasons of convenience. 

As has been already shown in [19], a particular choice of the realistic residual two-body 
interaction potential has almost no impact on the hnally calculated mean value and variance 
a of , with the overwhelming contribution to a coming from the choice of the single¬ 
particle basis size. Therefore, we perform the calculations here using G-matrix based only 
on the Bonn-CD nucleon-nucleon potential. 

It is well known that the residual interaction is an effective interaction depending on the 
size of the single-particle (s.p.) basis. Hence, when the basis is changed, the interaction 
should be modihed as well. For each nucleus in question three single-particle bases are 
chosen as described in Ref. [19] with the smallest set corresponding to Ihu particle-hole 
excitations, and the largest to about Ahoj excitations. The s.p. energies are calculated with 
the Coulomb corrected Woods-Saxon potential. 


A. Parameter adjustment 

In QRPA and RQRPA there are three important global parameters renormalizing the 
bare residual interaction. First, the pairing part of the interaction is multiplied by a factor 
Qpair whose magnitude is adjusted, for both protons and neutrons separately, such that the 
pairing gaps for the initial and hnal nuclei are correctly reproduced. This is a standard 
procedure and it is well-known that within the BCS method the strength of the pairing 
interaction depends on the size of the s.p. basis. 

Second, the particle-hole interaction block is renormalized by an overall strength param¬ 
eter Qph which is typically adjusted by requiring that the energy of the giant GT resonance 
is correctly reproduced. We find that the calculated energy of the giant GT state is almost 
independent of the size of the s.p. basis and is well reproduced with Qph ~ 1. Accordingly, 
we use Qph = 1 throughout, without adjustment. 

Third, an important strength parameter Qpp renormalizes the particle-particle interaction 
(the importance of the particle-particle interaction for the f3 strength was recognized first in 
Ref. [29], and for the f3f3 decay in [25]). The decay rate for both modes of f3f3 decay is well 
known to depend sensitively on the value of Qpp. This property has been used in [19] to £x 
the value of Qpp for each of the s.p. bases so that the known half-lives of the 2vj3j3 decay are 
correctly reproduced. 

Such an adjustment of Qpp, when applied to all multipoles has been shown in [19] to 
remove much of the sensitivity to the number of single-particle states, to the NN potential 
employed, and even to whether RQRPA or just simple QRPA methods are used. This is in 
contrast to typical conclusion made in the recent past [24, 30, 31] that the values of M'^'^ 
vary substantially depending on all of these things. 

We believe that the 2v decay rate is especially suitable for such an adjustment, in par¬ 
ticular because it involves the same initial and final states as the Oz/ decay. Moreover, the 
QRPA is a method designed to describe collective states as well as to obey various sum rules. 
Both double-beta decay amplitudes, Oz//?/? and 2uj3j3, receive contributions from many in¬ 
termediate states and using one of them for hxing parameters of QRPA seems preferable. 
We will elaborate this point in the next section. 

While the above arguments are plausible, the issue of adequacy of the QRPA method in 
general, and the chosen way of adjusting the parameters, should be also tested using suitable 
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simplified models that allow exact solution. One such test was performed recently [32], It 
involved two shells of varying separation, and a schematic interaction. From the point of 
view of the present work, the model suggested that QRPA is an excellent approximation of 
the exact solution. However, it turned out that the present method of parameter adjustment 
was not able to eliminate fully the effect of the higher, almost empty, shell on the Qu(3(3 decay 
matrix element. It is not clear whether this is a consequence of the schematic nature of the 
model, or of some more fundamental cause. In any case, the present results suggest that, 
within QRPA and RQRPA, the effects of the far away single particle states can be indeed 
eliminated, or at least substantially reduced. 

It is well known that the calculated Gamow-Teller strength is larger than the experimental 
one. To account for this, it is customary to ‘quench’ the calculated GT matrix elements. 
Formally, this could be conveniently accomplished by replacing the true value of the coupling 
constant gA = 1.25 by a qnenched valne 5 ^^ — 1-0. It is not clear whether similar phenomenon 
exists for other mnltipoles, besides J'^ = I"*". To see the dependence on the chosen valne, 
we nse in this work both the nnqnenched and qnenched valne of the axial cnrrent conpling 
constant gA = 1-25 and gA = 1-0, respectively (for all mnltipoles). The matrix elements M'^'^ 
calculated for the three s.p. bases and a fixed gA are relatively close to each other. As in [19], 
for each nncleus the corresponding average matrix elements (averaged over the three 

choices of the s.p. space) is evalnated, as well as its variance a. These qnantities (with the 
valne of a in parentheses) are shown in Table I, colnmns 4 and 5. Two lines for each nncleus 
represent the resnlts obtained with gA = 1-25 (the npper one) and gA = 1-0 (the lower one). 
One can see that not only is the variance snbstantially less than the average valne, bnt the 
results of QRPA, albeit slightly larger, are quite close to the RQRPA values. Furthermore, 
the ratio of the matrix elements calculated with different gA is closer to nnity (in most cases 
they differ only by ~20%) than the ratio of the respective gA sqnared (1.6 in onr case). The 
reason for snch a partial compensation of the gf^-dependence is that the experimental 
for gA = 1-0 is larger than for gA = 1-25 after adjusting in both calcnlation the gpp 
separately to the experimental 2z//5/5-decay transition probability. Gorrespondingly, one gets 
smaller adjusted value of gpp leading to larger calcnlated M^'^. Thns, with the adopted choice 
of parameter hxing the resnlting Oz//?/? decay rate depends on the adopted gA markedly less 
than the naive scaling g\ that would suggest a change by a factor of 2.44. 

Naturally, the 2vj3j3 half-lives are known only with some uncertainty. To see how the ex¬ 
perimental error in affects the calculated M'^'^, the derivatives at the experimental 
valne of are calculated. Finally, the errors indnced by the experimental nncertainties 
in Eexp. = given in colnmn 6 of Table I. 

Another nncertainty is related to the treatment of the short range nncleon correlation. 
The adopted form of the Jastrow-type factor /{ru), Eq. ( 8 ), is based on the work [28] where 
a range of the valnes of the parameter 71 is given. We vary that parameter 71 from 0.9 to 
1.2 (72 is fnlly determined by 71 ). The calculated dependence of M'^'^ is shown in Fig. 1 for 
^®Ge and ^°'^Mo. One can see that over this range of 71 valnes the calculated results differ by 
only abont 10 % from the ones corresponding to 71 = 1.1 nsed in the standard calculations. 

Since is not affected by the short range repnlsion. Fig. 1 simply shows that the QuPP 
matrix elements M'^'^ are not very sensitive to reasonable variations of the parameter 71 

[33]. 

Gombining the average (M' with the phase-space factors, the expected half-lives (for 
RQRPA and {m^) = 50 meV, the scale of nentrino masses snggested by oscillation experi¬ 
ments) are also shown in Table I (colnmn 7). 
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The information collected in Table I is presented in graphical form in Fig. 2. There the 
averaged nuclear matrix elements for both methods and both choices of qa are shown along 
with their full uncertainties (theoretical plus experimental). 

The entries for ^^°Te are slightly different from the corresponding results in Ref. [19]. This 
is so because in the present work we use the 2i>(3(3 lifetime tentatively determined in the 
recent experiment [34], while in [19] we used the somewhat longer 2vj3j3 lifetime based on 
the geochemical determination [35]. 


B. Preliminary discussion of the procedure 

One can qualitatively understand why our chosen procedure stabilizes the M'^'' matrix 
elements as follows: The matrix elements involve only the I"*" (virtual) states in the 
intermediate odd-odd nucleus. The nuclear interaction is such that the Gamow-Teller cor¬ 
relations (spin one, isospin zero pairs; after all the deuteron is bound and the di-neutron is 
not) are very near the corresponding phase-transition in the channel (corresponding to 
the collapse of the QRPA equations of motion). The contributions of the I"*" multipole for 
both modes of the (3(3 decay (2z/ and Ou) depends therefore very sensitively on the strength 
of the particle-particle force, parameterized by gpp. On the other hand, the matrix 

element, due to the presence of the neutrino propagator, depends on states of many multi¬ 
polarities in the virtual intermediate odd-odd nucleus. The other multipoles, other than I’*', 
correspond to small amplitudes of the collective motion; there is no instability for realistic 
values of gpp ~ 1.0. Hence, they are much less sensitive to the value of gpp. In the four panels 
of Fig. 3 we show the gpp dependence of the contribution of the I"*" multipole on one hand 
and of all the other multipoles added on the other hand. Remembering that the nominal 
value of gpp is ~ 1.0, one can clearly see the large difference in the corresponding slopes. 
By making sure that the contribution of the I"*" multipole is fixed, we therefore stabilize the 
value. The fact that RQRPA essentially removes the instability becomes then almost 
irrelevant thanks to the chosen adjustment of gpp. The effect of stabilization against the 
variation of the basis size is seen clearly if one plots directly M'^'^ versus Fig. 4. Also 
the obtained multipole decompositions of M'^'' plotted in Fig. 5 show the essential stability 
of the partial contributions against variation of the basis size. (Note that in the case of 
our approach gives rather small matrix element and its uncertainty is even compatible with 
zero value. The smallness of that matrix element is due to a large negative contribution of 
the I'*' multipole, see Fig. 5). 

Since the 2v lifetime depends on the square of there is an ambiguity in choosing 
the sign of the and hence the corresponding value of gpp. In this work, and in Ref. [19], 
always the solution corresponding to the smaller gpp value is chosen (positive with 
the present phase convention). It is worthwhile to justify such a choice. There are several 
reasons why the smaller gpp should be used. First, QRPA and RQRPA are methods designed 
to describe small amplitude excitations around the mean held minimum. Were we to choose 
the larger value of gpp, close or past the critical ‘collapse’ value, the method would be less 
likely to adequately describe the corresponding states. Second, as we will show later, by 
choosing the larger gpp the disagreement between the experimental and calculated rate of 
the single beta transitions from the lowest I"*" state in the intermediate nucleus would be far 
worse than with the smaller gpp. Moreover, as shown in Ref. [36], only with the smaller gpp 
can one successfully describe the systematic of single beta decay in a variety of nuclei. It 
follows from the study of Ref. [37] that choosing the larger value of gpp (i.e., the negative 
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sign of would lead to a complete disagreement with the systematics of single beta 

decays. Finally, there is also a pragmatic argument for such a choice. Only with it, the Oz/ 
nuclear matrix element becomes independent of the size of the single particle basis. Thus, 
this choice, admittedly ad hoc, removes the dependence on many more essentially arbitrary 
choices. 

One can see in Fig.5 that all multipoles (see Eq.(7)), with the exception of the and 
various very small entries, contribute with the same sign. This suggests that uncertainties 
in one or few of them will have relatively minor effect. It is instructive to see separately the 
effect of the sometimes neglected short-range repulsive nucleon-nucleon repulsion and of the 
induced weak nucleon currents. These effects are shown in Fig. 6 . One can see in that hgure 
that the conclusion of the relative role of different multipoles is affected by these terms. For 
example, in ^°'^Mo the 1 “ multipole is the strongest one when all effects are included, while 
the 2 “ becomes dominant when they are neglected. 

It is worthwhile to point out that one can display various contributions to M'^'^ in an 
alternative, perhaps more revealing, way. Instead of J'^ that represents virtual states in 
the intermediate odd-odd nucleus, one can decompose the result in terms of JT”, the angular 
momentum of the neutron-neutron (and simultaneously proton-proton) pair that undergoes 
the transition. In that case JT” = O"*" represents the ‘pairing only’ contribution, while JT” 7 ^ 0’*' 
come from the ground state correlations, i.e., contributions from higher seniority (broken 
pairs) states. In Fig. 7 we show that these higher seniority states contributions consistently 
have the tendency to cancel the ‘pairing only’ piece. Thus, the hnal is substantially 
less than the 77 = 0 ^ part only, signifying the importance of describing the ground state 
contributions properly. (This tendency is a well known effect, for both modes of the (3(3 
decay. It has been discussed, e.g. in Ref. [38].) 


III. UNCERTAINTIES OF THE 0u(3(3 DECAY MATRIX ELEMENTS 

One cannot expect that the QRPA-like and shell model calculations will lead to identical 
results due to substantial differences between both approaches. (However, we point out 
below that the results of the present approach and the shell model results differ relatively 
little whenever a comparison is possible.) Our goal in this section is to show that a majority 
of differences among various QRPA-like calculations can be understood. In addition, we 
discuss the progress in the held and possible convergence of the QRPA results. Based on 
our analysis, we suggest that it is not appropriate to treat all calculated 0z//5/5-decay matrix 
elements at the same level, as it is commonly done (see e.g., [18, 31]), and to estimate their 
uncertainty based on their spread. 

Here we shall discuss the differences among different published QRPA and RQRPA re¬ 
sults. We shall not consider the problem of the proton-neutron pairing in the double beta 
decay [39] and the self-consistent RQRPA (SRQRPA) calculations [40]. A systematic study 
of the 0z//5/3-decay matrix elements within the SRQRPA will be discussed in a forthcoming 
publication, where we show that within the considered approach of hxing nuclear structure 
input [19] the results agree well with those obtained within the QRPA and the RQRPA. 

Representative examples of the nuclear matrix elements calculated by different authors 
within the QRPA and RQRPA for nuclei of experimental interest are collected in Table 
H. In order to understand the differences between the entries, let us enumerate the main 
reasons leading to a spread of the published QRPA and RQRPA results: 



i) The quasiparticle mean field. Since the model space considered is finite, the pairing 
interactions have been adjusted to fit empirical pairing gaps based on the nuclear binding 
energies. This procedure is followed practically by all authors. However, it fails for closed 
proton (^^®Sn) and neutron (^^®Xe) shells. Some authors [31] modify the single particle 
energies (mostly spherical Woods-Saxon energies) in the vicinity of the Fermi surfaces to 
reproduce the low-energy quasiparticle spectra of the neighboring odd-mass nuclei. 

ii) Many-body approximations. The RQRPA goes beyond the QRPA by partially taking 
into account Pauli principle violation in evaluation of the bifermion commutators. Based 
on that, one might expect that the RQRPA is more accurate than the unrenormalized 
QRPA. Calculations within solvable models [41] support this conjecture. 

iii) Nucleon-nucleon interaction. In the 0z//5/5-decay calculations both schematic zero- 
range [42] and realistic interactions were considered. In Ref. [43] G-matrix of the Paris 
potential approximated by a sum of Yukawa terms was used. The interaction employed 
by the Tuebingen group has been the Brueckner G matrix that is a solution of the Bethe- 
Goldstone equation with Bonn (Bonn CD, Argonne, Nijmegen) one boson exchange po¬ 
tential. The results do not depend significantly on the choice of the NN interaction 
|19]. 

iv) The renormalization of the particle-hole interaction. This is achieved by scaling the 
particle-hole part of the (R)QRPA matrix by the parameter Pp^. That parameter is 
typically adjusted by requiring that the energy of some chosen collective state, often 
the giant Gamow-Teller (GT) resonance, is correctly reproduced. In Ref. [19] it was 
found that the GT state is almost independent of the size of the model space and is well 
reproduced with pph ~ 1. In Refs. [23, 39] Pph = 0.8 was chosen and in Ref. [44] pph = 1.3 
was considered for A > 114 nuclei. The sensitivity of results to the change of Pph were 
studied in Ref. [40] where it was shown that for 0.8 < Pph < 1.0 the change of the RQRPA 
0z//5/5-decay matrix element does not exceed 10%. 

v) The renormalization of the particle-particle interaction. This is achieved by scaling the 
particle-particle part of the (R)QRPA matrix by the parameter Ppp. When the early 
QRPA calculation were performed, only a limited information about the experimental 
2z//5/3-decay half-lives was available. Thus, in Ref. [42] a probable Ppp window was 
estimated from j3^ decays of semi-magic nuclei. In another work [43] systematics of the 
pp-force were investigated by an analysis of the single j3^/EC decays. Alternatively, 
in many works Ppp = 1.0 was chosen [23, 39, 40, 44, 45]. Nowadays, the 2z//5/5-decay 
has been observed in ten nuclides, including decays into two excited states [46]. Two 
recent papers used the 2z//5/5-decay half-lives to fix the strength of the particle-particle 
interaction of the nuclear Hamiltonian. In Ref. [30] it was used only for the J'^ = 
channel leaving the particle-particle strength unrenormalized (i.e., Ppp = 1.0) in other 
channels. In our previous paper [19], and in the current work, Ppp is adjusted so that 
the 2z//5/5-decay rate is correctly reproduced. The same Ppp is used for all multipole 
channels of the particle-particle interaction. Since information that can be used to adjust 
the Qpp value exists only for the 1^ channel, we cannot make a separate adjustment for 
other multipoles. Our economical choice is then to use the same Ppp for all multipoles. 
That makes the hamiltonian as simple as possible, and preserves the relative strength of 
different multipoles of the realistic starting point interaction. On the other hand, some 
authors prefer to fix the Ppp value to the f3~ decay transition of the ground state of the 
intermediate nucleus [31]. Such procedure can be carried out, however only for three 
nuclear systems (A=100, 116 and 128) where I’*' is the ground state of the intermediate 
nucleus. Since forbidden decays are less well understood than the allowed ones it is 
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difficult to rely on them in the case of other nuclear systems. (We discuss in more detail 
the differences between these two approaches of hxing the parameter Qpp in the next 
section.) 

vi) The size of the model space. A small model space comprising usually of two major shells 
was often used in the calculation of the nuclear matrix elements [30, 31, 42, 43, 44, 45]. 
Signihcantly larger model space has been used in Refs.[19, 23] consisting of hve major 
shells. As a rule the results obtained for the large model space are reduced when compared 
to those for the small model space for the same value of Qpp ~ 1. The dependence of the 
calculated nuclear matrix elements on the size of the model space has been studied in 
Ref. [30] and in [19]. Rather different conclusions using the same procedure of hxing Qpp 
was reached. The origin of such discrepancy is not understood. In Ref. [30] relatively 
small differences between model spaces were considered and a large effect on the Qu(3(3- 
decay matrix element was found. In particular, by adding the N = 2 shell to the N = 3—4 
oscillator shell model space the 0z//5/5-decay matrix element of decreased by a factor 
of about 3 (see Table II). Thus, in Ref. [30] the single-particle levels lying far from the 
Fermi surface seem to inhuence strongly the decay rate. In contrast, in Ref. [19] and in 
the present work results essentially independent on the size of the basis have been found 
even with signihcantly diherent sizes of the model space. 

vii) The closure approximation. The 0z//5/5-decay matrix elements were usually calculated 
using the closure approximation for intermediate nuclear states [30, 42, 43, 45]. Within 
this approximation energies of intermediate states {En — Ei) are replaced by an average 
value (< En — Ei >k. 10 MeV), and the sum over intermediate states is taken by closure, 

\Jn X ^nl — This simplihes the numerical calculation drastically. The calcula¬ 
tions with exact treatment of the energies of the intermediate nucleus were presented in 
Ref. [19, 23, 39]. The ehect of the closure approximation was studied in details in Ref. 
[47]. It was found that the diherences in nuclear matrix elements are within 10%. This is 
so because the virtual neutrino has an average momentum of ~ 100 MeV, considerably 
larger than the diherences in nuclear excitation energies. 

viii) The axial-vector coupling constant g^. The axial-vector conpling constant or in other 
words, the treatment of qnenching, is also a sonrce of diherences in the calculated nuclear 
matrix elements. The commonly adopted valnes are qa = 1-0 [42] and qa = 1-25 [19, 23, 
39, 43, 45]. However, as shown in Table II, if Qpp is hxed to the 2z//5/5-decay half-life, the 
ehect of qa modihcation is smaller, of order 10 %. 

ix) The two-nucleon short-range correlations (s.r.c.). In majority of calculations the short- 

range correlations between two nucleons are taken into account by mnltiplying the two- 
particle wave fnnctions by the correlation fnnction [28]: /(r) = 1 — (1 — 72 ^^) with 

7 i = 1.1 fm“^ and 72 = 0.68 fm“^. The sensitivity of resnlts to the change of this 
parameters have been discussed in the previous section. It is also known that the Qu(3(3- 
decay matrix element is more ahected (reduced) by the s.r.c. at higher Qpp values. We 
note that the s.r.c. were not taken into acconnt within the approach which was developed 
in Ref. [44] and used in the recent publication [31]. For the realistic Ppp values neglecting 
s.r.c. would lead to an increase in by factor of about 2 . 

x) The higher order terms of the nucleon current. The momentum dependent higher order 
terms of nucleon current, namely the induced pseudoscalar and weak magnetism, were 
hrst considered in connection with the light neutrino mass mechanism of the 0z//5/5-decay 
in [23]. Their importance is dne to the virtual character of the exchanged neutrino, with 
a large average momentnm of ~ 100 MeV. The corresponding transition operators have 
different radial dependence than the traditional ones. It is worth mentioning that with 
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modification of the nucleon current one gets a new contribution to neutrino mass mech¬ 
anism, namely the tensor contribution. The corrections due to the induced pseudoscalar 
nucleon current vary slightly from nucleus to nucleus, and in all cases represent about 
30% reduction. In the present calculation the effect of higher order terms of nucleon 
current is taken into account. 

xi) The finite size of the nucleon is taken into account via momentum dependence of the 
nucleon form-factors. Usually this effect on the Oufifi-decaj matrix elements associated 
with light neutrino exchange is neglected [42, 45], since it is expected to be small. In the 
presented calculation it is taken account by assuming phenomenological cutoff (see [23] 
and references therein). We found that by considering the hnite nucleon size the value of 
M'^'" is reduced by about 10%. In Refs. [31, 44] the nucleon current is treated in a different 
way than in other double beta decay studies, namely it is evaluated from the quark 
level using relativistic quark wave functions. However, the momentum dependence of 
corresponding nucleon form-factors is not shown, so a comparison is difficult. Comparing 
the present results with those of Ref. [44] for the same nuclear structure input suggests 
that the agreement might be achieved only if a very low cutoff is introduced. 

xii) The overlap factor of intermediate nuclear states. This factor is introduced since the two 
sets of intermediate nuclear states generated from initial and hnal ground states are not 
identical within the considered approximation scheme. A majority of calculations uses a 
simple overlap factor X'’X^ — Y'’Yfi inspired by the orthogonality condition of RPA states 
generated from the same nucleus. However, the double beta decay is a two-vacua problem. 
The derivation of the overlap matrix within the quasiboson approximation scheme was 
performed in Refs.[36, 48]. It was shown that the overlap factor of the initial and hnal 
BCS vacua is an integral part of the overlap factor of the intermediate nuclear states 
[36, 49]. This BCS overlap factor, about 0.8 for spherical nuclei, is commonly neglected 
in the calculation of nuclear matrix elements. If one hts Qpp to the experimental 2vj3j3- 
transition probability, the effect of the BCS overlap factor is signihcantly reduced. Hence 
it is neglected in the present work. This factor has been found to be important in the 
calculation within the deformed QRPA. 

xiii) The nuclear shape. Until now, in all QRPA-like calculations of the dufi-decay matrix 
elements the spherical symmetry was assumed as the majority of nuclei of experimental 
interest are nearly spherical. The effect of deformation on the 0z//5/5-decay matrix elements 
has not been studied as of now. Recently, the 2z//5/5-decay matrix elements were calculated 
within the deformed QRPA with schematic forces [36]. It was found that differences in 
deformation between initial and hnal nuclei have a large effect on the 2z//5/5-decay half-life. 
One could expect that a similar mechanism of suppression of nuclear matrix elements is 
present also in the case of the 0z//5/5-decay. It goes without saying that a further progress 
concerning this topics is highly desirable. The deformed QRPA might be the method 
of choice for the description of double beta decay of heavy nuclei like ^^°Nd, or ^®'^Gd. 
However, the deformation of many medium heavy nuclear systems are also noticeable (see 
Table H of [36]). The main difficulty of the deformed RQRPA calculation is the fact that 
by going from the spherical to the deformed nuclear shapes the number of conhgurations 
increases drastically. 

The present Table H, is organized as follows: i) The results with and without inclusion of 

higher order terms of nucleon current are separated. We point out that in Refs. [31, 44, 50] 

the contribution from the weak-magnet ism was taken into account. It is, however, negligible 

and we present their results without this contribution, ii) We indicate those matrix elements. 
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which differ signihcantly from the present ones or other previous calculations (denoted as 
SK-01 [30] and CS-03 [31]). We shall discuss this problem in detail below, hi) Within the 
sub-blocks the nuclear matrix elements are presented in ascending order following the year 
they appeared. 

Many nuclear matrix elements in Table II were calculated from the published ratios of 
Fermi and Gamow-Teller contributions and the absolute values of presented 

in units of fm“^, i.e., not scaled with nuclear radius [39, 44, 45, 50] We assumed that 
R = with ro = 1.1 fm. However, Ref. [43] used ro = 1.2 fm, so we rescaled the results 

in order to compare them better with ours. The differences between many calculations are 
understandable just from the way Qpp was hxed, the considered size of the model space, the 
inclusion of the short-range correlation, the way the hnite nucleon size is taken into account 
and other minor effects. Some nuclei, like ^'^°Mo, exhibit more sensitivity to these effects 
than others, such as '^®Ge or ®^Se. That is conhrmed by our numerical studies. 

The calculations of the 0z//5/5-decay matrix elements by Givitarese and Suhonen [31] 
(denoted as GS-03 in Table II) deserves more comments. The authors performed them 
within the approach suggested in Ref. [44], (SKF-91) which employs the nucleon current 
derived from the quark wave functions. In this approach the two nucleon short-range- 
correlations are not taken into account. But, the contribution from higher order terms of 
nucleon current is studied. Unlike the present results the authors of Ref. [44] claim that 
contribution from the induced pseudoscalar coupling is a minor one and they do not include 
it. At the same time, they find that the contribution from the weak-magnetism is negligible 
[44], in agreement with our result and with the l/rUp scaling. 

In Ref. [44] the Qpp was taken to be unity. The extension of this work for the case when 
Qpp is adjusted to reproduce the single /5-decay decay amplitudes was presented in Ref. [50] 
(AS-98) where also the effect of adjustment of the single particle energies to agree better 
with the spectroscopic data of odd-mass nuclei (AAaws —98) was studied. Since the adopted 
value of qa used for hxing Qpp was not given, in Table II we present the corresponding M'^'^ 
for both Qa = 1.0 and qa = 1-25. (From the related article [52] it seems qa = 1-0 was 
considered.) 

In Ref. [31] (GS-03) the nuclear matrix elements are calculated in the same way as in 
[50], however, the obtained results differ signihcantly (see Table II) from each other. For 
some nuclei the difference is as large as a factor of two. There is no discussion of this there 
or in the later Refs. [31, 52]. It is worth noticing that the largest matrix element in [31] 
is found for the 0z//5/5-decay of ^^®Xe. This disagrees with the results of other authors (see 
Table II). The reduction of the 0z//5/5-decay of the ^^®Xe is explained by the closed neutron 
shell for this nucleus. A sharper Fermi surface leads to a reduction of this transition. 

Altogether, the matrix elements of Refs. [31, 50] are noticeably larger than the present 
ones. Most of that difference can be attributed to the neglect of the short range nucleon- 
nucleon repulsion and of the higher order effects in the nucleon weak current in these papers. 

As the above list of 13 points shows, there are many reasons, some more important than 
others, that might cause a difference between various calculated 0z//3/5 matrix elements. 
Glearly, when some authors do not include effects that should be included (e.g. the short 
range correlations or the higher order terms in the nucleon current) their results should be 
either corrected or convincing arguments should be given why the chosen procedure was 
adopted. Other effects on the list are correlated, like the size of the model space and the 
renormalization of the particle-particle interaction. Again, if those correlations are not taken 
into account, erroneous conclusion might be drawn. Yet other effects are open to debate. 
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like the quenching of qa or the adopted method of adjusting gpp. In our previous work [19], 
and in the present one, we show that our chosen way of renormalization removes, or at least 
greatly reduces, the dependence of the final result on most of the effects enumerated above. 

Our choice of Qpp adjustment, to the experimental 2z//5/3 rate was used also in Refs. [30, 51]. 
Let us comment on their results. Note, that the higher order corrections (h.o.c.) to the weak 
nucleon current were neglected in those works, hence one expects a ~ 30% discrepancy right 
away. 

In Ref. [51] only the (3(3 decay of ^®Ge was considered within both the QRPA and RQRPA. 
Only one s.p. basis was used, corresponding to the small one in our notation. The result 
was = 4.59 (3.88) within QRPA (RQRPA). If one takes into account that h.o.c. as a 
rule reduce by about 30% and the fact that ~ 10% larger nuclear radius was used, one 
concludes that these matrix elements should be multiplied by ~ 0.6 in order to compare 
them with our calculations. This results in = 2.76 (2.34) which are in a very good 
agreement with our calculated = 2.68 (2.41) for this case. Hence QRPA and RQRPA 
gave quite similar results in that case. 

In contrast, the conclusions in Ref. [30] are quite different from ours; the authors found 
signihcant dependence of their results on both the nuclear model and on the s.p. basis size, 
in spite of the adjustment of Qpp. The reason for these differences is unknown. In an attempt 
to understand the origin of them we comment here on features that, in our judgment, require 
further discussion. 

1) We have already pointed out that the for '^®Ge in Ref. [30] changes with changing 
the single-particle basis within QRPA - from 4.45 (9 levels) to 1.71 (12 levels), or more than 
2.5 times. This was obtained by adding just the deep-lying sd shell, a surprising result. 

2) The last of the three papers contains most information and we can use their Figs. 1-14 

with and Figs. 15-25 with , in order to compare the numbers in the hrst (QRPA) 
and the last (SQRPA) columns of their Table 3. We were able to reproduce the entries in 
the QRPA column, but almost all entries in the SQRPA column, apart from those for ®^Se 
and seem inconsistent. The most striking difference is for ^^®Xe. By inspecting Fig. 

7 in [30], one hnds that the appropriate gpp for SQRPA is about 1.08 (lines a and b). Now, 
by going to Fig. 24 one can see that the corresponding valnes of are abont 2.5 (lines c 
and d) while the entries in Table 3 are 0.98 and 1.03. 

3) Some of the resnlts in [30] (see Fig.l, lines a and b, SQRPA; Fig.2, lines a and b, 
SQRPA, Fig.5, lines c and d, QRPA), Fig.7, lines c and d, QRPA) contradict the generally 
accepted conclnsion that for all RPA-like approaches the fnnctions M‘^'^{gpp) calculated with 
the larger basis cross zero faster than the ones obtained with the smaller one. Thns, given 
these apparent inconsistencies, it is difficult to draw any conclusions from the comparison 
of onr and Ref. [30] results. 

As pointed ont above, one cannot expect a perfect agreement between the present resnlt 
and the large scale shell model resnlts [54]. The approximations are different, and the shell 
model calculations do not inclnde several multipoles (which typically enhance the matrix 
elements, see Fig. 5). Yet, remarkably, the present results and the pnblished large scale shell 
model results agree with each other considerably better than the various entries in Table II. 
So, nsing the pnblished valnes, [55], we find M'^’^ valnes 1.5, 2.1, 1.1, and 0.7 for '^®Ge, ®^Se, 
^^'^Te and ^^®Xe, while RQRPA gives 2.4, 2.1, 1.5, and 0.7-1.0 for the same nnclei. However, 
it appears that the more recent shell model resnlts [56] for ^^°Te and ^^®Xe give larger values 
while at the same time overestimating to some degree the rate of the 2v(3(3 decay. In any 
case, we find this comparison enconraging. 
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IV. FURTHER DISCUSSION OF THE PARAMETER ADJUSTMENT 


Ideally, the chosen nuclear structure method should describe all, or at least very many, 
experimental data and do that without adjustments. As described above that is not the 
case of QRPA or RQRPA. The interaction used is an effective interaction, and various 
parameters {gpair, dph, dpp) are adjusted. In particular, the parameters gpair and gpp are 
adjusted on the case-by-case basis in the present approach. Even then the method is not 
able to describe well all relevant weak transitions. In particular, it is sometimes impossible to 
describe simultaneously the 2uPP decay rate as well as the (5~ and (3^/EC matrix elements 
connecting the I"*" ground states of the intermediate nucleus with the ground states of the 
final and initial nuclei (^'^°Mo is a well known example of this problem, see e.g. [57]). 

Empirically, the transitions through the ground state of ^°'^Tc, ^^®In, and ^^®I seemingly 
account for most of the 2vj3j3 matrix element (this is the so called ‘Single State Dominance’ 
[53]). Thus it appears that these single f3 transitions are particularly relevant. Based on such 
considerations, Suhonen [52] suggested that the f3~ matrix element is more suitable source 
for the gpp adjustment than the 2vj3j3 decay. Below we explain why we prefer the chosen 
method of parameter adjustment. On the other hand, we also argue that the adjustment 
based on the (3~ transitions might give results not far from ours, provided the independence 
on other adjustments can be proved. 

First, it is not really true that the hrst state is the only one responsible for the 2vj3j3 
decay. This is illustrated for the cases of ’^®Ge and ^°'^Mo in Fig. 8. Even though for ^'^°Mo 
the first state contributes substantially, higher lying states give non-negligible contribution. 
And in ^®Ge many I"*" states give comparable contribution. Thus, to give preference to the 
lowest state is not well justihed, the sum is actually what matters. At the same time, the 
dilemma that the (3~ and j3^/EC matrix elements move with gpp in opposite directions 
makes it difficult to choose one of them. It seems better to use the sum of the products of 
the amplitudes, i.e. the 2z//5/5 decay. 

At the same time, the contribution of the multipole to the Ouf]/] matrix element and 
the corresponding 2u/3/3 matrix element are correlated, even though they are not identical, 
as shown in Fig. 8. Making sure that the 2u/3/3 matrix element agrees with its experimental 
value constrains the part of the QuPP matrix element as well. 

In Fig. 9 we show the running sum contributions to the Oz//?/? matrix elements in ^®Ge and 
ioomo, separated into multipoles, and for the total. Such a sum, even for the component, 
is rather different that the similar staircase for the 2i'(3(3] there is no single state dominance. 
Thus, again, it is not obvious that it is best to choose any one particular state or transition 
for the adjustment. Moreover, while we have demonstrated that adjusting gpp to the 2z//5/5 
rate removes the dependence on other parameters, a similar proof was not given in [52]. 
Finally, as seen in Fig. 3 the slope of the multipole component of the Qvj3j3 matrix 
element far exceeds the slope of the other multipoles. Hence, making sure that the whole 
I"*" multipole contribution is correct is crucial. 

Thus, we prefer the method of adjustment used in the present work. That is so not only 
for the reasons shown above, but also because it is much more general as stressed already, 
while the ground state of the odd-odd nucleus has only in ^'^°Tc, ^^®In, and ^^®I. 

However, as demonstrated in the last hgure. Fig. 10, the adjustment proposed in Ref. 
[52], namely to choose the parameter gpp based on the f3~ decay of the intermediate I"*" state, 
would not be drastically different compared with the procedure used in the present work. 
As one could see, the resulting gpp are similar (but not identical) and the dependence on 
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the s.p. basis will be also reduced (That important feature was not demonstrated in Ref. 
[52], unfortunately). However, as the upper panel shows, it is difficult to describe the other 
component of that decay, the j3^/EC amplitude. That is an obvious drawback of the QRPA 
method; it is never meant to describe in detail properties of non-collective states. But that 
is less relevant for the description of integral quantities that depend on sums over many 
states. 

As pointed out earlier, we always choose the Qpp corresponding to the positive M'^’^ value. 
In the upper panel of Fig. 10 one can clearly see that were we to choose the other possibility, 
i.e., the Qpp corresponding to the negative the disagreement with the single beta decay 
would be considerably worse. It is important also to notice that if Pauli principle violation is 
restored (e.g. within the RQRPA) one hnds that the solution corresponding to the negative 
value of M‘^‘' is out of physical interval of gpp. This has been conhrmed also in a schematic 
model by presenting the solution of the QRPA with full inclusion of the Pauli exclusion 
principle [41]. As pointed out above, it was also found recently [36] that for negative value 
of M‘^'^ the correspondence with particle-particle strength from systematic studies of the 
single beta decay [37] is not achieved. 

In this section we have summarized the arguments why we believe that the procedure of 
adjustment used in the present work is preferable to the procedure advocated in Ref.[52]. 
At the same time, we suggest that there is no fundamental difference between the two; both 
are used to £x the fast varying contribution of the 1^ multipole. The bulk of the matrix 
element is associated with the other multipoles, and their effect is much less dependent on 
relatively small variations of the parameter Qpp. 

As we stressed already above, a substantial part of the differences in the calculated values 
of the Qu(3(3 decay matrix elements has its origin not in the choice of Qpp or other parameters, 
but in the neglect of the short range nucleon-nucleon repulsion and of the induced weak 
currents in some papers, and their inclusion in other papers (like the present one). Until a 
consensus on the treatment of these physics issues is reached, the differences in the calculated 
matrix elements cannot be avoided. 


V. SUMMARY AND CONCLUSIONS 

We have shown that the procedure suggested in our previous work. Ref. [19], is applicable 
to essentially all nuclei with known 2vj3j3 decay lifetimes . Adjusting the strength of the 
particle-particle neutron-proton force Qpp in such a way that the experimental 2z//5/5 decay 
rate is correctly reproduced removes much of the dependence on the size of the single-particle 
basis and whether QRPA or RQRPA is used. Here we also show that the quenching of the 
axial current matrix elements, parameterized by the reduction of the coupling constant qa-, 
also leaves the resulting 0z//5/5 matrix elements almost unchanged; they become insensitive 
to the variations of parameters describing the short-range nucleon-nucleon correlations as 
well. Thus, the resulting {)vj3j3 matrix elements acquire well dehned values, free of essentially 
arbitrary choices. We also present arguments while we believe that the chosen procedure of 
adjusting the interaction is preferable to other proposed ways of adjustment. 

We then summarize many published QRPA and RQRPA results and discuss their sim¬ 
ilarities and differences. We show that in most, albeit not all, cases these differences can 
be understood. We present an exhaustive list of reasons why individual calculated Oz//?/? 
nuclear matrix elements, evaluated within QRPA or RQRPA, might differ from each other. 
That list then can serve as a guide to readers of the past and future papers devoted to the 
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subject. Comparison between the results of different QRPA/RQRPA calculations would be 
facilitated if authors of future publications specify in detail what choices of explicit and 
implicit adjustable parameters they made, and discuss the dependence of their result on 
their particular choice. We believe that by following these suggestions a consensus among 
the practitioners of QRPA/RQRPA could be reached and most of the spread between the 
calculated nuclear matrix elements, that causes much confusion in the wider physics com¬ 
munity, would be shown to be essentially irrelevant. To reach a convergence of the results 
obtained using QRPA/RQRPA is clearly just an important step on the way to reliable and 
correct Qu(3(3 decay nuclear matrix elements. Exploring the structure of the intermediate 
odd-odd nuclei by the charge exchange reactions would create an opportunity to test the 
nuclear models more thoroughly. Further progress in nuclear shell model calculations and in 
the exploration of the exactly solvable models could also point the way towards the ultimate 
solution of this important problem. 


Acknowledgments 

Discussions with Jon Engel are gratefully acknowledged. The work of F. S. and V. R. 
was supported in part by the Deutsche Forschungsgemeinschaft (grants 436 SLK 17/298, 
TU 7/134-1 and FA67/28-2, respectively). We thank also the EU ILIAS project under the 
contract RII3-CT-2004-506222. The work of P. V. was supported by the U.S. DOE, Stanford 
University, SLAG and KIPAC. 


[1] Y. Fukuda et ah, Phys. Rev. Lett. 81, 1562 (1998); Y. Ashie et al, hep-ex/0501064. 

[2] B. T. Cleveland et ah, Astrophys. J. 496 , 505 (1998). 

[3] J. N. Abdurashitov et ah, Phys. Rev. C 60 , 055801(1999); W. Hampel et ah, Phys. Lett B 
447 , 127 (1999). 

[4] S. Fukuda et ah, Phys. Lett B 539 , 179 (2002). 

[5] Q. R. Ahmad et ah, Phys. Rev. Lett. 87 , 071301 (2001); Q. R. Ahmad et ah, Phys. Rev. Lett. 
89 , 011301(2002). 

[6] K. Eguchi et ah, Phys. Rev. Lett. 90 , 021802 (2003); T. Araki et ah, hep-ex/0406035. 

[7] R. D. McKeown and P. Vogel, Phys. Rep. 394 , 315 (2004). 

[8] M. C. Gonzales-Garcia and Y. Nir, Rev. Mod. Phys. 75, 345 (2003). 

[9] W. Grimus, hep-ph/0307149. 

[10] S.M. Bilenky, C. Giunti, J.A. Grifols, E. Masso, Phys. Rep. 379, 69 (2003). 

[11] P. Langacker, hep-ph/0411116. 

[12] A. Faessler and F. Simkovic, J. Phys. G 24, 2139 (1998). 

[13] J.D. Vergados, Phys. Rep. 361 , 1 (2002). 

[14] S. R. Elliott and P. Vogel, Annu.Rev.NuchPart.Sci. 52,115(2002); S.R. Elliott, Nuch Phys. 
Proc. Supph 138, 275 (2005). 

[15] S. R. Elliott and J. Engel, J. Phys. G 30 , R183 (2004). 

[16] S. J. Ereedman and B. Kayser, physics/0411216 (see also http://www.aps.org/neutrino/ for 
more details) 

[17] G. Aalseth et ah, hep-ph/0412300. 

[18] J.N. Bahcall, H. Murayama, and C. Pena-Garay, Phys. Rev. D 70, 033012 (2004). 


16 



[19] V.A. Rodin, A. Faessler, F. Simkovic and P. Vogel, Phys. Rev. C 68, 044302 (2003). 

[20] Oi^/3/3 decay can be also mediated by the exchange of various heavy particles. For a review of 
this possibility see [12]. Relation to the lepton flavor violation has been suggested in [21]. 

[21] V. Cirigliano, A. Kurylov, M. J. Ramsey-Musolf, and P. Vogel, Phys. Rev. Lett. 93 , 
231802(2004). 

[22] S. Pascoli and S. T. Petcov, Phys. Lett. B 544, 239 (2002). 

[23] F. Simkovic, G. Pantis, J.D. Vergados, and A. Faessler, Phys. Rev. C 60 , 055502 (1999). 

[24] F. Simkovic, J. Schwieger, G. Pantis, and Amand Faessler, Found. Phys. 27 1275 (1997). 

[25] P. Vogel and M. R. Zirnbauer, Phys. Rev. Lett. 57 , 3148 (1986). 

[26] O. Givitarese, A. Faessler, and T. Tomoda, Phys. Lett. B 194, 11 (1987). 

[27] K. Muto and H. V. Klapdor, Phys. Lett. B 201, 420 (1988). 

[28] G. A. Miller and J. E. Spencer, Ann. Phys. 100, 562 (1976). 

[29] D. Cha, Phys. Rev. C 27, 2269(1983). 

[30] S. Stoica and H. V. Klapdor-Kleingrothaus, Eur. Phys. J. A 9 , 345 (2000); Phys. Rev. C 63 , 
064304 (2001); Nucl. Phys. A 694 , 269 (2001). 

[31] O. Givitarese and J. Suhonen, Nucl. Phys. A 729 , 867 (2003). 

[32] J. Engel and P. Vogel, Phys. Rev. G 69 , 034304 (2004). 

[33] Similar conclusions were reached by J. Engel earlier in a study involving an exactly solvable 
model. We appreciate his input. 

[34] G. Arnaboldi et al. Phys. Lett. B 557 , 167 (2003). 

[35] T. Bernatowicz et al, Phys. Rev. Lett. 69, 2341 (1992); T. Kirsten, H. Richter, and E. 
Jessberger, Phys. Rev. Lett. 50, 474 (1983). 

[36] F. Simkovic, L. Pacearescu, and A. Faessler, Nucl. Phys. A 733 , 321 (2004); R. Alvarez- 
Rodriguez, P. Sarriguren, E. Moya de Guerra, L. Pacearescu, A. Faessler and F. Simkovic, 
Phys. Rev. G 70 , 064309 (2004). 

[37] H. Homma et al.., Phys. Rev. G 54, 2972 (1996). 

[38] J. Engel, P. Vogel, Xiangdong Ji, and S. Pittel, Phys. Lett. B 225 , 5 (1989). 

[39] G. Pantis, F. Simkovic, J.D. Vergados, and A. Faessler, Phys. Rev. C 53 , 695 (1996). 

[40] A. Bobyk, W.A. Kaminski, and F. Simkovic, Phys. Rev. C 63 , 051301 (2001). 

[41] F. Simkovic, A.A. Raduta, M. Veselsky, and A. Faessler, Phys. Rev. G 61 , 044319 (2000). 

[42] J. Engel, P. Vogel and M. R. Zirnbauer, Phys. Rev. C 37 , 731 (1988). 

[43] K. Muto, E. Bender, and H.V. Klapdor, Z. Phys. A 334, 187 (1989). 

[44] J. Suhonen, S.B. Khadkikhar, and A. Faessler, Nucl. Phys. A 535, 509 (1991). 

[45] T. Tomoda, Rept. Prog. Phys. 54 , 53 (1991). 

[46] V.I. Tretyak and Yu.G. Zdesenko, At. Dat. Nucl. Dat. Tabl. 80, 83 (2002). 

[47] K. Muto, Nucl. Phys. A 577, 415c (1994). 

[48] F. Simkovic, G. Pantis and A. Faessler, Prog. Part. Nucl. Phys. 40, 285 (1998); Phys. Atom. 
Nucl. 61, 1218 (1998). 

[49] K. Grotz and H. Klapdor, Phys. Lett. B 157, 242 (1985). 

[50] M. Aunola and J. Suhonen, Nucl. Phys. A 643 , 207 (1998); Nucl. Phys. A 602 , 133 (1996). 

[51] K. Muto, Phys. Lett. B 391, 243 (1997). 

[52] J. Suhonen, Phys. Lett. B 607 , 87 (2005). 

[53] A. Garcia et ah, Phys. Rev. C 47 , 2910 (1993); H. Ejiri and H. Toki, J. Phys. Soc. Jpn. 65 , 7 
(1996); M. Bhattacharaya et ah, Phys. Rev. G 58 , 1247 (1998); O. Givitarese and J. Suhonen, 
Nucl. Phys. A 653 , 321 (1999); F. Simkovic, P. Domin, and S.V. Semenov, J. Phys. G 27 , 
2233 (2001); P. Domin, S.G. Kovalenko, F. Simkovic, and S.V. Semenov, Nucl. Phys. A 753 , 


17 



337 (2005). 

[54] E. Caurier, F. Nowacki, A. Poves, and J. Retamosa, Phys. Rev. Lett. 77, 1954 (1996). 

[55] E. Caurier, F. Nowacki, A. Poves, and J. Retamosa, Nucl. Phys. A 654, 973c (1999); E. Cau¬ 
rier, G. Martinez-Pinedo, F. Nowacki, A. Poves and A. P. Zuker, nucl-th/0402046. 

[56] F. Nowacki, private communication. 

[57] A. Griffiths and P. Vogel, Phys. Rev. G 46, 181 (1992). 

[58] H. Homma, E. Bender, M. Hirsch, K. Muto, H.V. Klapdor-Kleingrothaus, T. Oda, Phys. Rev. 
C 40 (1989) 540. 


18 



TABLE I: Averaged Ouf3f3 nuclear matrix elements and their variance a (in parentheses) 

evaluated in the RQRPA and QRPA. In column 6 the variance Eexp. of the 0j^/?/3-decay matrix 
element due to uncertainties in the measured 2z//3/3-decay half-live is given. and qa 

denote the 2j^/3/3-decay nuclear matrix element deduced from and axial-vector coupling 

constant, respectively. In column 7 the QvfSfd half-lives evaluated with the RQRPA average nuclear 
matrix element and for assumed {mpp) = 50 meV are shown. For ^^®Xe there are four entries; the 
upper two use the upper limit of the 2v matrix element while the lower two use the ultimate limit, 
vanishing 2v matrix element. ^^®Nd is included for illustration. It is treated as a spherical nucleus; 
deformation will undoubtedly modify its matrix element. 


Nuclear 


QA 


(M' 


^exp. 

to 

((m/ 3 / 3 ) = 

: 50 meV) 

transition 


[MeP-i] 

RQRPA 

QRPA 



[yrs] 




1.25 

0.15 ±0.006 

2.40(0.07) 

2.68(0.06) 

±0.05 


^•^^-0.16 

1027 



1.00 

0.23 ±0.01 

2.30 (0.04) 

2.48 (0.05) 

±0.05 


2 40+0-15 

1027 



1.25 

0.10 ±0.009 

2.12 (0.10) 

2.36 (0.09) 

±0.04 


'3-'30’-0.62 

1026 



1.00 

0.16 ±0.008 

1.91 (0.05) 

2.10 (0.07) 

±0.05 


0 1 0+0.64 

1026 

mzr 


1.25 

nil -1-0.03 
'-'•-‘--‘■-0.06 

0.31 (0.08) 

0.04 (0.10) 

- 1 - 0.20 

- 0.43 


1 

-*-•^^-0.96 

1028 



1.00 

'0.06 

0.43 (0.11) 

0.40 (0.02) 

-1-0.25 

-0.45 


D 77®^ 

' ^-0.48 

1028 

- 


1.25 

0.22 ±0.01 

1.16(0.11) 

1.28(0.09) 

±0.02 


1 4(3+0-30 
l.OU-o 23 

1027 



1.00 

0.34 ±0.015 

1.12 (0.09) 

1.24 (0.08) 

±0.02 


1.45l0:il 

1027 

- 


1.25 

0.12 ±0.006 

1.43 (0.08) 

1.56 (0.10) 

±0.03 


8-20lJ:°^ 

1026 



1.00 

0.19 ±0.009 

1.22 (0.07) 

1.31 (0.08) 

±0.02 


11.3^}:^ 

1026 

i2%Te - 

> 128^6 

1.25 

0.034 ± 0.012 

1.60 (0.11) 

1.73 (0.13) 

±0.09 


1 op:+0.38 
1.00-0,29 

1028 



1.00 

0.053 ±0.02 

1.37 (0.07) 

1.47 (0.05) 

±0.1 


9 p;9+0.52 
^•^^-0.40 

1028 

i30Te - 


1.25 

0.036lO;[]39 

1.47 (0.15) 

1.55(0.17) 

- 1 - 0.3 

-0.09 


0-871°;?^ 

1027 



1.00 

0.056lHi 

1.28 (0.08) 

1.36 (0.10) 

-1-0.27 

-0.08 


1 1 c:+0.23 
-‘-•-‘-^-0.38 

1027 

i36Xe - 


1.25 

0.030 

0.98(0.09) 

1.03(0.08) 



-*-•^4-0.30 

1027 



1.00 

0.045 

0.90 (0.07) 

0.94 (0.05) 



9 1 q+0-25 

1027 



1.25 

0 

0.73(0.09) 

0.77(0.10) 



0 09 + 1 .00 
^•^^-0.69 

1027 



1.00 

0 

0.63 (0.07) 

0.67 (0.08) 




1027 

1^0 Nd - 


1.25 

n n'7+0.009 
^•^^-0.03 

2.05 (0.13) 

2.25 (0.16) 

-1-0.07 

- 0.20 


0.92lO;26 

1026 



1.00 

n 11 -I-0.014 
'-’•-‘--‘--0.05 

1.79 (0.05) 

1.96 (0.04) 

-1-0.06 

-0.19 


1 91 H“0.32 
-‘■*^■'-- 0.10 

1026 
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TABLE II: Neutrinoless double beta decay matrix element M'^’^ calculated within QRPA and 
RQRPA approaches. The way of fixing the particle-particle interaction strength is indicated: a 
given value, to 2jv/3/3-decay half-life, to f3 decays. The nuclear radius is R = Notation: 

EVZ-88 = Engel, Vogel, Zirnbauer[42] (for a'l = -390), MBK-89 = Muto, Bender, Klapdor[43], 
T-91 = Tomoda [45], SKE-91 = Suhonen, Khadkikhar, Eaessler[44], PSVE-96 = Pantis, Simkovic, 
Vergados, Paessler[39], AS-98 = Aunola, Suhonen[50], results for Woods-Saxon (WS) and adjusted 
Woods-Saxon (AWS) bases, SPVE-99 = Simkovic, Pantis, Vergados, Paessler[23], SK-01 = Stoica, 
Klapdor-Kleingrothaus[30] (the second line by SK-01 results corresponds to slightly larger model 
space), CS-03 = Civitarese, Suhonen[31]. 

Ref. Method vq gpp qa 

[fm] '^^Ge 82Se 

Without higher order terms of nucleon current 









Differences understandable 



EVZ-88 

QRPA 

1.1 

/3 

1.0 

1.3 

1.0 

0.8 

1.8 


2.4 

2.2 

1.0 


MBK-89 

QRPA 

1.1 

/3 

1.25 

3.84 

3.59 


1.94 


3.93 

3.18 

1.45 

5.57 

T-91 

QRPA 

1.1 

1.0 

1.0 

2.86 

2.59 


3.14 


2.55 

2.22 

1.27 

3.47 





1.25 

3.97 

3.60 


4.30 


3.53 

3.07 

1.74 

4.80 

SKE-91 

QRPA 

1.1 

1.0 

1.00 

3.37 

2.72 




3.17 

2.94 

1.71 






1.25 

4.55 

3.71 




4.24 

3.95 

2.31 


PSVE-96 

QRPA 

1.1 

1.0 

1.25 

3.04 

2.23 

2.41 

1.09 

0.94 

2.48 

2.33 

1.55 


SPVE-99 

RQRPA 

1.1 

1.0 

1.00 

4.05 

3.82 

2.24 

4.58 

2.86 

3.38 

2.87 

1.20 

5.15 





1.25 

3.60 

3.40 

1.99 

4.12 

2.58 

2.96 

2.50 

1.02 

4.51 

present 

QRPA 

1.1 

2uPP 1.25 

3.35 

2.95 


1.83 


2.32 

1.98 

1.30 

3.10 









Results discussed 




SK-01 

QRPA 

? 

2uPP 

? 

4.45 

5.60 

4.16 

5.37 

3.99 

4.84 

4.73 

1.69 







1.71 

4.71 

2.75 

3.81 

2.85 

3.43 

3.77 

1.35 



RQRPA 

? 

2v[5[5 

? 

3.74 

4.30 

3.01 

4.36 

3.61 

4.29 

4.55 

1.57 







1.87 

2.70 

2.72 

3.40 

3.39 

2.83 

3.00 

1.02 


AS'vv5-98 

QRPA 

1.1 

/3 

1.00 

3.98 

3.69 

2.88 


2.21 


4.62 

2.49 






1.25 

5.30 

4.93 

3.85 


2.93 


6.15 

3.34 



QRPA 

1.1 

/3 

1.00 

4.85 

3.61 

3.70 


3.97 


3.81 

2.15 






1.25 

6.44 

4.82 

3.96 


5.25 


5.05 

2.84 


CS-03 

QRPA 

? 

/3 

1.25 

3.33 

3.44 

3.55 

2.97 

3.75 


3.49 

4.64 








With higher order terms of nucleon current 


SPVE-99 

RQRPA 

1.1 

1.0 

1.25 

2.80 

2.64 

1.49 

3.21 

2.05 

2.17 

1.80 

0.66 

3.33 

present 

QRPA 

1.1 

2v[5[5 1.00 

2.48 

2.10 

0.40 

1.24 

1.31 

1.47 

1.36 

0.94 

1.96 





1.25 

2.68 

2.36 

0.04 

1.28 

1.56 

1.73 

1.55 

1.03 

2.25 


RQRPA 

1.1 

2ul3(5 1.00 

2.30 

1.91 

0.43 

1.12 

1.22 

1.37 

1.28 

0.90 

1.79 





1.25 

2.40 

2.12 

0.31 

1.16 

1.43 

1.60 

1.47 

0.98 

2.05 
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FIG. 1: Nuclear matrix element M'^'^ as function of parameter 71 of the func 
used to account for the two-nucleon short range correlations. Calculation p' 
with Qnn = 0.99 C^^Ge), 1.21 and the small single particle level scher 
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FIG. 2: Average nuclear matrix elements and their variance (including the error coming 

from the experimental uncertainty in for both methods and for all considered nuclei. For 

^^®Xe the error bars encompass the whole interval related to the unknown rate of the 2vj3j3 decay. 
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0.6 0.7 0.8 0.9 1.0 1.1 1.2 0.7 0.8 0.9 1.0 1.1 1.2 
^PP ^PP 

FIG. 3; Contribution of the 1"*“ and all the other (summed) multipoles to the matrix elements M'^'^ 
as a function of the parameter gpp. The results were obtained within the RQRPA approach for the 
large model space. The dots indicate the gpp values that reproduce the known 2u decay rate. 








FIG. 4; Dependence of the matrix elements M'^'^ on M‘^'^ originating from the variation of the 
parameter gpp only. The results were obtained within the QRPA approach. The vertical lines show 
the experimental M'^’^ values used for the gpp adjustment. 
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FIG. 6: The effects of higher-order terms of nucleon currents (h.o.c.) and of the nucleon-nucleon 
short range repulsion (s.r.c.) on the multipole distribution of the Oz^/3/3-decay matrix element in 
the QRPA. The left panels were calculated with 9(13) single-particle levels and the right panels 
with 21 single-particle levels. 
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FIG. 7: Running sum of the 0u(5(5 decay matrix element as a function of the angular momentum 
77 of the nn and pp pair that undergoes transition. 
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FIG. 9: Running sum of the 0i^/3/3-decay matrix element for and and their multipole 

contributions as a function of excitation energy E^x = En — (Ei + Ef)l2. Calculations were 
performed within the QRPA (small model space). 
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FIG. 10: The upper panel: The matrix element associated with the f3~ transition from the initial 
nucleus (A,Z) to the I"'' ground state of the intermediate nucleus (A,Z-|-1) as function of Qpp for 
and The middle panel: The same as upper panel for the f3~ transition from 

the ground state of (A,Z-|-1) to ground state of final nucleus (A,Z-|-2). The lower panel: The 
2v(5fd-dec&y matrix element of and as a function of gpp. The experimental 

values are indicated by thin horizontal lines {qa = 1-25 is considered). 
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